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Abstract 

A formalism for determining energy eigenstates of cylindrical lead salt quantum wires in the 
multiple-band envelope-function approximation is developed. Electron energy dispersion for quan- 
tum wire subbands within the conduction and valence bands is found. 
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The effect of quantum confinement on electron and liole states in splierical quantum 
dots and cylindrical quantum wires of III-V and II- VI semiconductor compounds having a 
complex valence-band structure was described within the multiband envelope-function ap- 

n 

proximation more than 20 years ago [Jj. Description of electronic structure of spherical 
IV- VI semiconductor quantum dots using similar approximation followed soon 2|. Re- 
cently nanostructures based on IV- VI semiconductor compounds such as lead salts received 
much attention due to their potential for applications in solar cells and as infra-red de- 
tectors. Quasi-zero-dimensional nanostructures of lead salts have become the subject of 
wave-function engineering [s], 0], and there appeared an interest in quasi-one-dimensional 
nanostructures of IV- VI semiconductors Q] . However, one is forced to accept the fact that 
a recent attempt [5|] to describe electronic structure of a cylindrical lead salt quantum wire 
within a multiband envelope-function formalism is not quite correct. This Brief Report seeks 
to compensate for this deficiency. 

The conduction and valence band extrema in lead salt semiconductors (PbSe, PbS) occur 
at the L-points of the Brillouin zone. Electron spectrum near the L-point taking into account 
only the two closely lying conduction and valence bands and neglecting band anisotropy can 
be described by the spherical Dimmock model [gI-SI. In this model the electron wave 
function is written as 

m = u\L~,)+v\Li), (1) 

where |Lg ) and \Lq) describe the Bloch functions while u{t) and 'O(r) are the spinors slowly 
varying with coordinates and satisfying the equations jol 

(f-«,A) -iP{(tV 

Here cr/3 (/3 = x, y, z) are the Pauli matrices, ac, a^, Eg, and P are parameters of the model 
and E is the electron energy. 

We will first construct linearly independent solutions of Eqs. (|2]) having cylindrical sym- 
metry and describing electronic states in a bulk semiconductor. These solutions can be 
characterized by the energy, E, momentum, kz along the quantum wire axis, and projec- 
tion, M of the total angular momentum on the wire axis. Then we will impose the boundary 
condition of the four-component envelope wave function vanishing on the cylindrical surface 
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of the quantum wire. In the cyhndrical coordinates {p, (p, z) this boundary condition is to 
be imposed at p = R, where R is the radius of the cyhndrical quantum wire. 

The high symmetry of the problem allows us to separate variables and pin down a set of 
good quantum numbers {E, kz, and M). After this Eqs. ([2]) can be considered as a system 
of four coupled ordinary linear differential equations of second order with respect to the 
variable p. As we are only interested in solutions of Eqs. ([2]) finite at p = 0, there are at 
most four such linearly independent solutions of Eqs. ([2]). 

Let us look for a solution of Eqs. ([2]) in the form 

'Ce'(^^-'/'^^JM-i/2{kp)- 

De^^''^^l')^JM+i/2{kp), 
where Jn(x) is the Bessel function of order n, k is the wave number of transverse motion, 

and A, B, C, and D are the coefficients to be determined. Substitution of Eqs. ([3]), (jl]) into 

Eqs. ([2]) leads to a system of the following algebraic equations: 

PkzA-iPkB-{a^k'^ + a^kl + E + Eg/2) C = 0, (5) 

iPkA-PkzB-{a^k'^ + a^kl + E + Eg/2) D = 0, (6) 

(ttc k'^ + ackl-E + Eg/2) A + P k^C - i P k D = , (7) 

(a, k'^ + ack1-E + Eg/2) B + i P k C - P k^ D = . (8) 
The condition that this system of algebraic equations has a non-trivial solution yields 

k'' + kl = E + A, (9) 
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where 



[E{ay - a,) - Egjay + a,)/2 - P^Y + a,ay{4E^ - Eg) 

2acay 

Equations (|5]), (|6]) allow one to express the coefficients C and D in terms of A and B: 

P{kzA-ikB) 
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ayk^ + aykl + E + Eg/2 
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P{ikA-k^B) 



a^k'^ + a^kl + E + Eg/2 ' 
We see that a natural choice of the two hnearly independent solutions will be to set either 
A 01 B equal to zero. The resulting solutions (which turn out to be orthogonal) take the 
form 
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a^k'^ + a^kl + E + Eg/2 
The two remaining solutions of Eqs. ([2]) can be sought in the form 



(10) 

(11) 
(12) 

(13) 



um{p,^,z) = e 
vm{p,V,z) = e 
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■Fe^^"^-^l^)^lM-i/2M 

Ge^^^^'+y^)^lM+i/2M 
where In{x) is the modified Bessel function of order n while C, D, F, and G represent a new 

set of the coefficients to be determined. Substitution of Eqs. (JH]), (fT5|) into Eqs. ([2]) yields 



P KC - I P K D + {a^ - a^kl - E - Eg/2) F = 0, 



(16) 



- i P kC - P k, D + {a^ - a^kl - E - Eg/2) G = , (17) 
(oc kl- a^K^ - E + Eg/2) C + P k, F - i P kG = , (18) 
{ac kl- acK^ - E + Eg/2) D - i P k F - P k,G = . (19) 
The condition that this system of algebraic equations has a non-trivial solution yields 



(20) 



Equations f lT6|) . f|T7|l allow one to express the coefficients F and G in terms of C and D: 

^ _ P{-k,C + iKD) 

— otv kl — E — Eg/2 ' 
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FIG. 1: (Color online) Energy of electron states in PbSe cylindrical quantum wire at A;^ = as a 
function of the wire radius. 
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By setting zero either C or D we arrive to the following solutions of Eqs. (|2]) 



Wm(p,^,^) = Ce 



kz 2 



ay — kl — E — Eg/2 
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k"^ — a^kl — E — Eg/2 
The four solutions we constructed are not mutually orthogonal but they are linearly 
independent. If one requires that their linear combination (with the coefficients A, B, 
C, and D) vanishes at p = R then one will obtain a system of four homogeneous algebraic 
equations on these coefficients. The condition that this system has a non-trivial solution will 
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FIG. 2: (Color online) Energy dispersion of electron subbands within the conduction band of a 
PbSe cylindrical quantum wire of radius i? = 20 A. By dashed lines are shown subbands with the 



quantum number ric = 2. 



lead to the dispersion equation determining the allowed energy values of electrons confined 
in a cylindrical quantum wire. The positive (negative) values of energy describe the states 
in the conduction (valence) band. The dispersion equation takes the form 



a^k^{k + K ) JM-i/2{kR) JM+i/2{kR) Im-i/2{k.R) hi+i/2{.i^R) 

+ k{a^K^ - a^kl~ E - Eg/ 2) JM-i/2{kR) Im+i/2{i^R) 
+K (a^ k'^ + a^kl + E + Eg/2) JM+i/2{kR) /M-i/2(fi:i?)] (25) 
k {a^ - a^kl - E - Eg/ 2) JM+i/2{kR) Im-i/2{k.R) 
-K {a^ k'^ + a^kl + E + Eg/2) JM-i/2{kR) Im+i/2{i^R) 
In order to better understand the structure of this equation let us first consider the case 



oi kz = Q. In this limit each of the spinors wy, v 



(0 .-,(0 
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1,2,3,4) has only one non-zero 



component. Therefore, for every i the corresponding bispinor solution has two non-zero 
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FIG. 3: (Color online) Energy dispersion of electron subbands within the valence band of a PbSe 
cylindrical quantum wire of radius = 20 A. By dashed lines are shown subbands with the 
quantum number = 2. 

components. Taking into account the parity of the Bloch functions {L"^), the solutions with 
i = 1 and i = 3 have in the hmit kz = the parity (—1)*''^+^/^, while the solutions with 
i = 2 and i = 4 have the parity (—1)^^^^/^. The first term of Eq. fl25|) vanishes in the limit 
kz = 0, and Eq. (!25|) reduces to a product of the two square brackets. The first (second) 
square bracket is responsible for the solutions with the parity (—1)*'^^^/^. This situation is 
analogous to the case of spherical quantum dots Is]. 

When ^2 7^ then each bispinor solution we constructed has three non-zero components. 
In this case the first term in Eq. (!25|) does not vanish and factorization of the dispersion 
equation becomes impossible. This observation is in contradiction with the conclusions of 



Ref. 
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Because of the Kramers degeneracy Eq. f l25|) must be invariant under the change of M 
to —M. This invariance is guaranteed by the properties of the Bessel functions: J-n{x) = 
(-1)" J„(x); lJx) = lJx). 
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Let us illustrate our results by numerical calculations. In Fig. 1 are shown energies of 
electron states in a PbSe cylindrical quantum wire at fc^ = as a function of the wire radius. 
Only dependences for the five lowest subbands in the conduction band and five uppermost 
subbands in the valence band are shown. All the material parameters are taken from Ref. 
The electron states in the conduction (valence) band are characterized by the projection, 
M of the total angular momentum onto the wire axis and by the main quantum number Uc 

In Fig. 2 (Fig. 3) is shown energy dispersion of electron subbands within the conduction 
(valence) band of a PbSe cylindrical quantum wire of radius i? = 20 A. The solid lines refer 
to the subbands with the main quantum number nc(v) = 1 while the dashed lines refer to the 
subbands with the main quantum number nc{v) = 2. The dependences shown in Figs. 2, 3 
are nearly parabolic. 

To summarize, we have studied electronic structure of cylindrical lead salt quantum wires 
within the multi-band envelope-function formalism. We have found the energy dispersion 
for subbands within the conduction and valence bands of the quantum wire. When applied 
to real structures, the parameters of the model may be modified for each particular L-valley 
to reflect its orientation with respect to the quantum wire growth direction and taking into 
account the valley anisotropy. 
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